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Abstract 

We present an extension of the Edwards model for conformations 
of individual chain molecules in solvents in terms of fractional Brown- 
ian motion, and discuss the excluded volume effect on the end-to-end 
length of such trajectories or molecules. 

1 Introduction 

Individual chain polymers in good solvents are typically modelled by trajecto- 
ries of random walks, or - in the continuum limit - by Brownian paths. Such 
models by themselves however do not take into account that self-crossings 
of these paths should be suppressed, this "the excluded volume" effect will 
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make the trajectories less curly and more extended. Fractional Brownian 
paths have been suggested as a heuristic model for such swelling, or on the 
other hand for polymers in a collapsed state [2], but a more proper model 
would be based on self-avoiding random walks, or on a weight factor which 
penalizes self-crossings, such as in the continuum Edwards [3] [6] [20] [21] 
[2"2"] [2"5] or the discrete Domb- Joyce [5] model. 

The ensuing swelling of the molecular conformations is given by the Flory 
index [8] [9] which describes the scaling of the end-to-end distance as a func- 
tion the number of monomers. It has been extensively studied both in the 
(chemical) physics and the mathematics community. The physics literature is 
characterized by structural intuition and far-reaching predictions, the math- 
ematical results are less far-reaching but provide the high reliability charac- 
teristic of the mathematical approach. Both are too vast to be quoted here, 
we refer for this to recent reviews [12] |18j . 

In the present paper, after a few words on fractional Brownian motion 
fBm, we shall see that one can extend to fBm the Edwards model of Brownian 
paths with exponentially suppressed self-intersections, a mathematical exis- 
tence proof has been established recently [UJ. In the third part of the paper 
we generalize some by now classical arguments from the physics literature to 
explore what the Flory index might be in the fBm case. 

2 The fBm Edwards Model 
2.1 Fractional Brownian Motion 

Fractional Brownian motion on M. d , d > 1, with "Hurst parameter" H e (0, 1) 
is a d- dimensional centered Gaussian process B H = {B H (t) : t > 0} with 
covariance function 

E(B»(t)Bf(s)) = 5 f(t 2H + s 2H -\t-s\ 2H ), i,j = l,...,d, s,t>0. 

For H = 1/2 it is ordinary rf-dimensional Brownian motion B. We refer 
to the recent monographs by Biagini et al. [T] and by Y. Mishura [16]; for 
self- intersection local times of fBm see Hu and Nualart [T5] . 
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2.2 The Edwards Model 

Self-repelling Brownian paths for a time interval < t < I can be modelled 
via a "Gibbs factor" to suppress self-intersections: 

G = I exp f-g J^ds j dtS {B{s) - B(t)) \ . 

Technically one defines this expression as a limit, using 
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in particular 



Z = lim E exp [ -g ds dt5 £ (B(s) - B(t)) 



o jo 



if this quantity is well defined; otherwise a renormalization is required, as, 
more generally, in Theorem 2.2 below. 

Recently, generalizing an argument of Varadhan [20J, this was extended 
in pi] to 

G = iexp (-9 J ds f Q dt5 (B H ( s ) ~ B H (t))^ > 

as follows. 

Theorem 1 The Edwards model is well defined for all H < l/d, with 

G = iexp (-9 J o ds jT dtS (B H (s) - B H (t))^j . 
Theorem 2 For H = l/d and g sufficiently small 

G = lim i- exp jf jf dt5 e - ) j , 



with 



Z £ = E ^exp ^-0 jf ds jf cM £ - 5 H (t) )^ 



zs well-defined. 
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3 The Flory Index 



When the number N of monomers of a polymer becomes large one expects 
its end-to-end length R to scale [10] 



R{N) ~ N v . 

For (fractional) Brownian motion the root-mean-square length 

R= (B h (N) 2 ) 

is scaling with 

v = H. 

But the excluded volume effect makes the paths and polymers swell: the 
end-to-end length increases. For the Brownian motion case there is the fa- 
mous Flory formula 

based originally on a mean field argument. Since its proposal by Flory [7] 
[9], numerous methods were invoked to put it on a more solid mathematical 
basis, a process which has up to now been fully successful in the case d — 1 

To obtain what may be considered as a first guess of a similar formula for 
fBm we shall return to the modest beginnings, generalizing Fisher's original 
argument [7J [8] (see e.g. the review given in McKenzie pSJ) to the case at 
hand. 



3.1 The Fisher Argument 

A partition function Z(R) for a freely jointed chain of N segments for which 
the end-to-end length has fixed modulus R is given by 

Z(R) = aR d - 1 exp( ' lU ~ 



2N" 

and leads to a free energy 

F 1 = -lnZ~ w -(d-l)lnR. 
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Instead of such a chain a continuous model is that of a Brownian trajectory 
from time zero to time N, for which one computes 

E (6 (B (N) - R)) = (2nN)- d/2 exp (-^j . (1) 
For the fBm case this formula generalizes to 

E (s (B" (N) -R)) = (2«N™)-" 2 exp (-^/) (2) 
from which we see that N — > N 2H , and hence we should consider 

Z{R) = aR d ~ 1 e W {- Wm ) 

i.e. 

dR 2 

F 1 ^-\nZ^^ m -(d-l)\nR. 
For the repulsive excluded volume energy of fBm paths x with x(N) = R, 

F2 = — In E x ^_g f exp f — # J J dt6 (x(s) — x(t))\ 

dimensional considerations and mean field arguments [15] suggest 

N 2 

F 2 ~ const.— 

Maximizing 

F(N,R)=F 1 (N,R) + F 2 (N,R) 

with regard to R leads to 

<i-R d — 1 . t0 j 1 

= — - — - ctm sLN*ir+-\ 

Assuming that the 2nd term is negligible one finds 

R d + 2 ~ N 2H+2 

i.e. 

R ~ N Vh 

with 

fA\ 2H + 2 (*\ 
v H {d) = -j^. (3) 



Remark 3 A polymer model with B(t 2H ) instead of B H (t) would produce 
the same expression as in (TJ|) ; hence also the same Flory index, but would 
not share the homogeneity implied by the stationary increments of B H (-). 

3.2 The Critical Dimension 

The derivation of vjj is evidently heuristic and needs validation. For this it 
is worth noting that for Brownian motion there is a critical dimension d c = 4 
defined by the fact that for d> d c there is no excluded volume effect, so that 
R scales like the unperturbed Brownian motion: 



We can ask for which dimension, more generally, the fBm Flory index will 
show no excluded volume effect from self-crossings, i.e. 

v H (d c ) = H. 

Inserting our ansatz (J3]) one finds 

Hd c = 2. 

and indeed it is known (Theorem 1.1 of Talagrand [19]) that <i-dimensional 
fBm has no double points iff 

Hd > 2, 

in other words, our vh predicts d c correctly. 

Remark 4 As a consequence, any Flory formula should be considered only 
up to the critical dimension, i.e. as long as there are double points and 
an excluded volume effect. Similarly, any prediction of % > 1 would be 
unphysical: the end-to-end distance cannot grow faster than the number N of 
monomers. In the case at hand this suggests for the one- dimensional case 



R~ N 



1/2 



i.e. 



v l/2 (A) = 1/2. 
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Note that for small H the scaling exponent as predicted would be strictly less 
than one while for the Brownian motion case 

ui/ 2 (l) = 1 
has been proven [Wj f23\/ . The infimum 

lim v H (l) = 2/3 

happens to be the scaling exponent of the myopic random walk fT^j . 

Remark 5 In the attached figure 1 the two red lines correspond to vn(d) = 1 
and to the critical dimension as a function of the Hurst index H , respectively. 
Above these the Flory index is unphysical. On the green lines vh is validated. 
(The existence proof of the fBm Edwards model in fTDJ works below the dashed 
line.) 

Remark 6 For fixed dimension d, any extension F(H) of the Flory formula 
to general Hurst indices H will have to obey 

F(~) = -A- (4) 
K 2 J d + 2 v ; 

for the usual Brownian motion (Flory-Fisher), and 

2 2 

at the critical point (Talagrand). Note that our ansatz (Gj) 

2H + 2 



F(H) = v H (d) 



d + 2 

is just the unique linear interpolation between those two values. 



3.3 A Recursion Formula 

For H = 1/2 Kosmas and Freed [H] derive a recursion formula 



(Here and in the following we label formulas from - or analogous to those in 
- the paper [14] by their numbers in that article, in double brackets.) 

The derivation of this formula is specific to the Brownian motion case 
and does not hold for general %. Hence in what follows we shall generalize 



their arguments which led to ( (4.13) ) to first obtain a valid recursion formula 
and then check whether it is satisfied by vg as given in 02]). 
We begin by considering 

Z (g, N) = E ^exp (-g ds jf dt5 (B H (s) - B H (t)]j^j . 

From the defining relation 

E (B H (s)B H (t)) = l - [s 2H +t 2H -\s- t| 2H ) 

we see that for a > the processes {B H (t) : t > 0} and [a~~ H B H (at) : t > 0} 
obey the same law. Making this substitution and a change of integration 
variables as = a, at = r we obtain 

(/ paN paN \ \ 

exp I -ga Hd - 2 / ds dt5 (B H (s) - B H (t)) J J 

= Z(a Hd ~ 2 g,aN). ((2.14)) 

Likewise we find for the mean-square end-to-end distance 

(R 2 ) ee [{B H (N))\ W (- g J N da£dtS (B H (s) - B H (*)))) 

(6) 

(/ paN paN 

(B H (aN)) 2 ex V l-ga Hd ~ 2 / ds dtS (B H (s) - B H {t)) 

= a~ 2H f (a Hd - 2 g,aN) ((2.18)) 
= N 2H f(N 2 ~ Hd g,l). ((2.19)) 

(Note the critical dimension Hd = 2 where (R 2 ) 1 ^ 2 ~ N H .) For large N 
one expects a power law behavior for the unknown function /, i.e. 

(R 2 ) ~ iV 2H (iV 2 - H ^) x ((2.22)) 
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with an exponent x to be determined. 

As a next step we restrict one coordinate of the positions, Xi(t) = Bf (t) 
to the interval [0, D] by inserting 



into (Ml. One obtains 



1 if Bf (t) G [0, D] for all t 
otherwise 



(R 2 ) D = Zd (l[ ,D]( B ?) {B H (N))\xp (-gj ds j dt6 (B H (s) - B H \t)) 

-2H ( of r aN r aN 

= ZD{g n) E l [0 , a n D] {Bf) (B H (aN)) 2 exp -a Hd ~ 2 g J da J dr5 (B H (a) - B 

= a~ 2H F(a H D, a Hd ~ 2 g, aN) = N 2H F(N~ H D, N 2 ~ Hd g, 1). ((4.3)) 

Now assume that asymptotically there is a dimensionless correction factor 
h for 



D 



It should grow as D becomes small which suggests a power law behavior for 
the function h: 

^° » ^ ( TPy) " ((4 - 5)) 

with y to be determined. As D approaches a minimal value Dq - approxi- 
mately the extension of a monomer (" Kuhn length" ) - the polymer becomes 
effectively (d — l)-dimensional: 

<i&i> « <iS> D « ^ <^> 1+ * ■ ((4-6)) 
This provides a relation between the end-to-end length for dimensions d and 
d — 1. To obtain from this a recursion relation, recall equation ((2.22) ): 

and introduce instead of x the (unknown) 

2v H (d) = 2H + x{2 - Hd) 
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so that 

(R 2 d ) = c d N 2v »^g^^r- ((4.8)) 

and 

, 2v H {d-l)-2H 

(RU) = c^N^-^g . ((4.9)) 



On the other hand from ( (4.6) ) we have 

{RU) « (R 2 d ) 1+ > = const.N 2 ^ 1 ^) g 2 ^^^) . (7) 
Comparing exponents in these two expressions we find 

v H (d-l) = VH (d)(l + ^j 

v H {d- l)-E v H (d) - E / y \ 
2-E-(d-l) 2- Ed V 2/ 

The first of these equations gives 

y = v H (d-l) 
2 v H (d) ' 

with this the second one becomes 

1 v H (d-l)-H 1 v H {d)-H 
v H {d-l)2 - H ■ (d-l) ~ v H (d) 2 -Ed 

i.e. this expression does not depend on the dimension d so that all the vn(d) 
are given in terms of e.g. vh(1) : 

1 v H (d)-H 1 v H (l)-E 



v H {d) 2- Ed v H {l) 2-H ' 
^ " (d -l)v H (l) + 2 -dH- (8) 



Proposition 7 

^ 2H + 2 

satisfies this recursion equation, with 

1 i>^(d) - # 1 



2 -Ed 2E + 2' 
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Remark 8 The standard Flory index (TJj) for H = 1/2 obeys the recursion 
formula. 



Remark 9 The recursion formula (TJj] implies the correct critical behavior, 
i.e. any solution will obey vn(d) = H for d = 2/H = d c , whatever the choice 
o/i>i/(l). To see this explicitly, insert d = 2/H and find 

v (-) = gzjM! = H 

Hy H> {2/H-l)v H {l)+0 

Remark 10 7/t>#(l) turned out to be equal to one for all H, the recursion 
formula would suggest 

an expression which then also produces the standard Flory formula for H = 
1/2, as well as the critical dimension d = 2/H. 



4 Summary 

The Edwards type model for self-repelling fBm now at hand will raise the 
question of how the end-to-end length of trajectories scales as a function of 
time (or "number of monomers"). The original Fisher argument, while crit- 
icized regarding its assumptions jTT] [4], provides a simple heuristic "deriva- 
tion" of the Flory formula which allows an extension to fBm. The obtained 
scaling law needs further verification; we note that it correctly predicts the 
critical dimension for which the excluded volume become negligible and obeys 
a recursion formula based on dimension reduction. The latter would provide 
a useful constraint on any alternate scaling laws. 



* * * 
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Figure 1: The domain of the Flory index 
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